Effect of normal current on kinematic vortices 
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Within the framework of time-dependent Ginzburg-Landau theory, we discuss an effect of the 
non-magnetic interaction between the normal current and the supercurrent in the phase-slip regime. 
The correction due to the current-current interaction is shown to have a transient character so 
that it contributes only as a system evolves. Numerical analyses for thin layers with no magnetic 
feedback show that the largest contribution of the current-current interaction appears near sample 
edges, where kinematic vortices reach majcimum velocity. 



I. INTRODUCTION 

Although time-dependent Ginzburg-Landau (TDGL) 
theory is justified only for slowly evolving systems, it 
provides qualitatively correct description of such phe- 
nomena as the ultrafast propagation of magnetic flux 
dcndritesiiSi, fast kinematic vortices^i^ or accelerated vor- 
tex motion involving vortex-antivortex annihilation^. Of 
course, in these fast processes the quasiparticles cannot 
achieve a local equilibrium distribution, meaning that 
non-equilibrium corrections to the TDGL theory become 
important. This was demonstrated by Vodolazov and 
Peeters^ who found a large deformation of the gap pro- 
file at the phase-slip centre in the case of a slow relaxation 
of quasiparticles. 

Vodolazov and Peeters have assumed isotropic distri- 
bution of quasiparticles (valid in the dirty limit) and care- 
fully treated the energy distribution using two coupled ki- 
netic equations for longitudinal and transverse branches. 
Their approach applies for a finite gap, the time deriva- 
tive of which acts as a force driving quasiparticles out of 
equilibrium. 

Here we shall discuss a complementary correction 
(valid in the pure limit) which takes into account a 
direction-dependent perturbation of the momentum dis- 
tribution of quasiparticles; such a perturbation appears 
as the normal current is created. The correction to the 
TDGL equation is found to be proportional to the scalar 
product of the normal current and the supercurrent <i 



indirectly via the magnetic field. The absence of any di- 
rect interaction between these two currents in the TDGL 
theory is not disturbing, because it is in agreement with 
an intuitive picture based on the two-fiuid model of a su- 
perconductor: taking the condensate as an independent 
fluid one expects it not to interact with the underlying 
crystal including its normal electrons. 

The absence of interaction between normal current 
and supercurrent is also supported by microscopic the- 
ories in the dirty limit4i^— These approaches, however, 
cannot be used to discuss the current-current interac- 
tion. To obtain practical equations, authors employ the 
isotropic approximatioiiSiSiii in some cases assuming in 
addition local equilibrium^. The isotropic distribution 
corresponds to the zero normal current, therefore any 
effect of the normal current on a formation of supercon- 
ducting gap is lost in this approximation. 

The effect of the normal current on the gap has been 
derived in ?| from the Thouless criterioni^ adapted to 
non-equilibrium Green functions. In [isj it was applied 
to the FIR conductivity of the Abrikosov vortex lattice 
and shown to explain a decrease of the real part of con- 
ductivity below the critical temperature, experimentally 
observed with ultrafast spectroscopy^'', while the TDGL 
theory predicts a small increase. Since the microscopic 
derivation is lengthy and technically demanding, in the 
appendix we provide a simple derivation of the inter- 
action of the condensate with the normal current using 
purely phenomenological arguments. 



A. Normal current in a superconductor 

Superconductors with freely-moving Abrikosov vor- 
tices or propagating dendrites have a finite resistivity; 
an electric field E' thus develops in them as the current 
is driven through. This electric field generates a normal 
current 

Jw = CTatE', (1) 

which adds to the supercurrent J5. In the TDGL theory 
the normal current and the supercurrent interact only 



B. Plan of paper 

The paper is organized as follows. In section |lT] we in- 
troduce the floating-kernel approximation, which is the 
TDGL theory extended by the interaction between the 
normal current and the supercurrent. In section Hill we 
show that this correction is of transient nature, being 
zero in any steady regime. To this end in section llll Al we 
perform a gauge transformation to express conveniently 
the interaction of the normal current and supercurrent in 
terms of time derivatives of the vector and scalar poten- 
tials. Consequently, this correction can be described in 
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terms of effective magnetic and electric fields as shown 
in section IIII Bl In section IIVI we apply the theory to 
the phase-slip regime in a strip made of thin supercon- 
ducting layers with negligible magnetic feedback. After 
rescaling so all quantities are dimensionless, in section lTVl 
we present results of a numerical simulation to demon- 
strate how the current-current interaction influences fast 
kinematic vortices in the phase-slip regime. Section |V] 
contains concluding discussion. In the appendix we indi- 
cate why the TDGL theory violates the longitudinal /- 
sum rule and show that addressing this problem with an 
intuitive two-fluid correction leads to the floating-kernel 
approximation. 



II. FLOATING-KERNEL APPROXIMATION 

Here we write down a closed set of equations forming 
the theory we term the floating-kernel (FK) approxima- 
tion. This becomes identical to TDGL theory in the coor- 
dinate system floating with normal electrons, because the 
normal current vanishes in this reference frame. When 
the normal current is accelerated, the condensate experi- 
ences inertial forces absent in the laboratory TDGL the- 
ory. 



A. Order parameter 

In the presence of normal current the time-evolution 
of the order parameter is described byiii^ 
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The right hand side has been derived by Kramer 
and Watts- Tobin three decades ago^^ Terms aip and 
/? 1-01^ V' are a standard part of Ginzburg-Landau theory. 
The kinetic energy in the left hand side has been ob- 
tained only recently!^ It applies to pure superconductors, 
when the Cooperon mass equals twice the electron mass, 
m* = 2m. In the appendix this normal-current correc- 
tion is deduced from the longitudinal /-sum rule. The 
TDGL equation^ obtains setting Jjv = 0. 

Away from the critical line the phase and amplitude re- 
lax with different rates. The parameter controlling this 
difference is C = 2TinAo/(?i|V'oP), where Aq and tAo are 
values of the BCS gap and GL function at given tem- 
perature in the absence of currents. Since in pure super- 
conductors TinAecs ^ ^, the correction C^j-f/'P can be 
large under realistic conditions. We use this relaxation 
rate in the numerical example. Close to the critical line 
this correction vanishes and one can use a simpler theory 
corresponding to the limit C -^^ of equation 

Our major concern will be the contribution of the nor- 
mal current, the Jjv-term. In ^ the kinetic energy de- 



pends on the difference between the velocity of the con- 
densate 
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where -0 = j-^le'-*^, and the mean velocity of normal elec- 
trons 
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The first term of ([2]) is thus the kinetic energy which 
must be yielded by a pair of normal electrons in order to 
join the condensate, in the reference frame floating with 
normal electrons. To distinguish the theory based on the 
velocity difference from the standard TDGL theory, we 
call equation ^ the fioating-kernel approximation^^ 



B. Two-fluid picture of current 

The derivative of the kinetic energy with respect to 
vector potential A defines the current operator. The 
correction to the normal current thus also appears in the 
supercurrent 
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with ns = IV'P being the density of Cooper pairs or the 
condensate density. That this supercurrent depends on 
the relative velocity of the condensate with respect to the 
normal background is desirable. According to Ohm's law 
([1]) all electrons move with the normal velocity v^. If a 
superconducting fraction moves with a different velocity 
vg, we must add the difference. 

One may quit the picture of relative motion and re- 
arrange the total current in the spirit of the two-fiuid 
model 



e*ns (vs - vjv) + envAr 
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In this rearrangement the supercurrent has the conden- 
sate velocity vs . The correction term has become a part 
of the normal current, where it reduces the density of 
electrons to the fraction of normal electrons. We have 
used that Cooperon charge is twice the electron charge 
e* = 2e. 

The necessity to reduce the normal current to the nor- 
mal fraction follows from the longitudinal /-sum rule. In 
the appendix we show that in order to achieve a con- 
sistent theory formulated in terms of a free energy, the 
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reduced normal current must be accompanied by changes 
in the free energy which lead to the floating-kernel ap- 
proximation. 



C. Scalar and vector potential 

The vector and the scalar potentials A and (f) yield the 
electric field 
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In some applications one should keep in mind that is 
thought of as a local electrochemical potential, and not 
the electrostatic potential. The vector E' is thus the 
driving force per electron rather than the Maxwell elec- 
tric field. Following the notation of Josephson we write 
E' rather than E as a reminder of this distinction. As is 
usual in the theory of superconductivity we call E' the 
electric field for brevity. 

Although the system has non-zero scalar potential, de- 
viations from charge neutrality arc so small that one can 
neglect them using the continuity equation in its station- 
ary form V ■ J = 0. Substituting for the normal current 
from ([!]), one finds the usual condition for the potential 



(8) 



We have used V • A = and assumed a homogeneous 
sample; Vctat = 0. To evaluate the vector potential we 
need the Maxwell equation 



V^A = -/io (Js + Jtv) 



(9) 



which is also in the stationary approximation to be con- 
sistent with the continuity equation. The set of equations 
(PHI), (0, and dZlini) is closed. 



III. TRANSIENT NATURE OF THE 
INTERACTION OF THE NORMAL CURRENT 
WITH THE SUPERCURRENT 

An overlap of the normal current and supercurrent ap- 
pears at the conversion layer at the junction of the super- 
conductor to a normal lead. Similarly, there is such an 
overlap at phase-slip centres in superconducting wires or 
at phase-slip lines in films. In this section we show that 
the normal current is purely transient and contributes 
only if the electric and magnetic fields change in time. 



A. Effective vector and scalar potentials 

The normal current enters the floating-kernel approxi- 
mation in two ways: in the kinetic energy of the equation 
(12) and in the supercurrent ([3]). In both cases Jjv and 



A appear together so that both are accounted for by a 
vector field 
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It is advantageous to describe the vector and scalar 
potentials in a symmetric form. To this end we express 
the normal current ([TJ via potentials 
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By the substitution ip — c ^'^ '^'^1^ -0 for a homogeneous 
sample, Vr = 0, we transform the GL equation ([2]) as 
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The supercurrent ([5]) then reads 



Js = -^Re ( -im - —Aeff ) i> 
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and other equations of the TDGL theory need not be 
written again since they depend only on the amplitude 

i^p = iv;|2. 

It is now clear that the system behaves as if the nor- 
mal electrons are driven by potentials (f> and A while the 
superconducting electrons experience effective potentials 
(j)cs and Aoff. 



B. Effective electric and magnetic fields 

The above effective potentials give rise to effective 
magnetic and electric fields. The transverse component 
of the normal current acts on the condensate via an ef- 
fective magnetic field 



dB 

Bcff = V X Aeff = B - T— . 
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(18) 
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The time variation of the normal current acts on the 
condensate via an effective electric field 

1 d c)E' 
E'cff = — — Aoff- V0cff = E'-r— . (19) 
c ot at 

In both effective fields the correction term vanishes 
in the stationary limit. The corrections following from 
the floating-kernel picture might thus become important 
in transient regimes or in systems driven by oscillating 
fields. The AC response of the Abrikosov vortex lattice 
has been discussed in !l3'|. Here we focus on vortices 
driven by a steady supercurrent. 

IV. FAST VORTICES 

When driven by large transport DC^ vortices move fast 
and the condensate is strongly reduced behind the vortex 
core since it needs a finite time to recover. Similarly, the 
condensate is stronger on the front side since the con- 
densate needs a finite time to dissolve or move away. In 
consequence, vortices are so deformed that their interac- 
tion becomes anisotropic leading to transition from the 
triangular Abrikosov lattice to rows perpendicular to the 
current. These rows are known as phase-slip lines. Along 
the phase-slip lines vortices move with velocities exceed- 
ing typical velocities of normal vortices by one to two 
orders of magnitude.'*'^ 

In this regime one can expect the floating-kernel cor- 
rections to play an important role for two reasons. First, 
in the vicinity of fast vortices time derivatives of the mag- 
netic and electric fields are large so that the effective 
magnetic field ([TS]) or electric field (fTOl) appreciably dif- 
fers from the true field. Second, across the phase slip line 
there is a large normal current since vortices are densely 
packed there and the condensate is suppressed. 

Above a critical current the resistivity of a supercon- 
ductor changes in dramatic steps corresponding to rear- 
rangement of vortices inside. Presently, this highly non- 
linear response of the system can be studied only numer- 
ically. We thus demonstrate the effect of floating-kernel 
corrections using numerical simulation. 

A. Dimensionless equations for thin films 

In films much thinner than the London penetration 
depth, the current has a negligible feedback effect on the 
magnetic field. We thus assume a constant homogeneous 
magnetic field perpendicular to the superconducting film. 
According to (|T8| the effective magnetic field coincides 
with the true magnetic field. Bog = B. It is thus con- 
venient to take a stationary vector potential, dA/dt = 
in which case from (|16p follows that Acff = A. We thus 
focus on the scalar potential and floating corrections to 
it. 

The vector potential is constant and the Maxwell equa- 
tion ^ is unused. We must solve the TDGL equation 



(|14p which in dimensionless units reads 
(V'-iA')'V'+(l-|V'f) ^' 

where the order parameter ip' = ip/ipo is scaled with the 
GL value ipQ ~ —a/ (3. The amplitude of relaxation rate 
thus reads 7 = Cipo- Dimension x' = x/^ scales with 
the GL coherence length — —fi^/2m*a. The vector 
potential scales with inverse distance; A' = {e*£,/ch)A. 
The time t' = tu/rQh is scaled with the GL time tgl = 
—T/a = Trh/8kB(Tc — T) and phase relaxation rate u 
specified below from the scaling of the equation for the 
scalar potential. The scalar effective potential 4>'^^ = 
(f'cs/va scales with the inverse time ipo — hu/e*TGi^. 
The supercurrent (jl7p simplifies to 

J's = Re[ip' {-iV ~ A')i;'] (21) 

with = (m*^/e*fi^Q) J5. Equation ([8]) for the scalar 
potential rescales as 

V^<j>' = V' -J's, (22) 

where cf)' is in the same units as (p'^g. We have used u = 
nh/AkBTcT with r from Eq. to make Eq. free of 
numerical factors. The temperature dependence of tql 
cancels with the condensate density i/'o ~ n{T — Tc)/Tc. 
The effective potential is 

0;« = ^'-2(l-|)A^'. (23) 

In the numerical study we use values w = 6, 7 = 10 and 
T — 0.75 Tc. The sample has size 20^ x 40f and contacts 
of length 4^ are centered at shorter sides. Our numer- 
ical code is a modified version of that used by Berdiy- 
orov, Milosevic and Peelers. Principally, the modifica- 
tion has consisted of the addition of the normal-current 
term. The reader interested in details related to bound- 
ary conditions can find them in their paper [3] or in a 
closely related p^ . 

B. Numerical results 

Figures [T] and [5] show the order parameter in a super- 
conducting film in a weak magnetic field H = Q.l Hc2 
driven by a current from normal contacts. Near the cen- 
ter of the strip the current is strong enough to enforce a 
rearrangement of vortices into closely packed lines, while 
near the end points the vortices remain isolated. The 
strip thus sports both kinetic vortices in lines and iso- 
lated Abrikosov vortices. 

The left-hand panels show dynamics evaluated within 
the TDGL theory, the central panels show the dynamics 
obtained from the fioating-kernel approximation. The 
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FIG. 1. Effect of the ffoating-kernel correction on tlie con- 
densate density in a 2-D sample: Low density is indicated 
by blue and maximum by red. The left-hand panel shows 
the TDGL theory while the central panel shows the floating- 
kernel approximation. The sample is in a weak magnetic field 
which breaks left-right symmetry. The DC is driven into it 
from normal-metal contacts creating visible conversion lay- 
ers at horizontal edges. Near the center of the sample there 
are two phase-slip lines with kinetic vortices, closer to the 
contacts one can find isolated Abrikosov vortices. The right- 
hand panel presents the floating- kernel correction rdcfi/dt. It 
is small in magnitude being a few percent of the scalar po- 
tential (j). There is no visible contribution due to conversion 
layers, as this would correspond to the transient nature of the 
floating-kernel correction. The largest contributions (strong 
red and blue) appear near the ends of the phase-slip lines, 
where kinetic vortices achieve the highest velocity.-^ 




FIG. 2. As in figure [T] but for a current increased by 15%: 
Again, there are no traces of the conversion layers in the 
floating-kernel correction. The only appreciable contribution 
results from the kinetic vortex just at the sample edge. In the 
upper phase-slip line one can see a four-leafed clover show- 
ing the quadrupole character of the floating-kernel correction 
Tdcj>/dt of a single kinematic vortex. Similar features can also 
be identified for other kinematic vortices inside the sample. 

overall behaviour in both approximations is very simi- 
lar. The floating-kernel correction to the scalar potential 
leads to slightly different motion of vortices as one can 
see from their positions after the same time of evolution 
has elapsed from initial conditions of zero current and 
zero magnetic field. 



The right-hand panel in figures [T] and [2] show the 
floating-kernel correction to the scalar potential. It is 
viforthy of note that there are no visible contributions near 
contacts, where the normal current is converted into the 
supercurrent and these two currents have a strong over- 
lap giving seemingly good conditions for their interaction. 
This demonstrates the transient nature of floating-kernel 
corrections, as argued in section Hill 

It should be noted that for the above regime the true 
scalar potential has dominant contributions from conver- 
sion regions and barely-visible features from phase-slip 
lines. The amplitude of the floating-kernel correction is 
only few percent on this scale. 

There are no visible contributions to the effective po- 
tential from isolated Abrikosov vortices. This is due to 
their low velocities^ which renders their contribution in- 
significant on the present scale. 

Near sample edges the floating-kernel correction to the 
scalar potential is dominated by kinematic vortices. This 
corresponds to the acceleration of kinematic vortices to 
velocities by one order of magnitude larger near edges 
than in the interior of the samplcj^!^ 

In spite of a weak signal from vortices deep inside the 
sample, one can see in the right-hand panel of figure [2] 
that each vortex brings a quadrupole floating-kernel cor- 
rection to the potential. This can be understood from the 
Bardeen-Stephen picture.-- They have shown that a vor- 
tex moving in the horizontal direction creates an electric 
field which drives the normal current through its core in 
the vertical direction. The corresponding scalar poten- 
tial is thus a dipole with vertical orientation. The time 
derivative of this potential due its horizontal translation 
has quadrupole symmetry. 




FIG. 3. Voltage across the sample as a function of time: The 
TDGL theory (black symbols) and floating-kernel approxi- 
mation (red line) yield similar mean voltage, VroGLC/Voi = 
0.0416 and VFK^/v'oi = 0.0440. In the floating-kernel ap- 
proximation the voltage oscillates faster than in the TDGL 
case. The voltage is a difference of potentials at centres of 
contacts. 

Figure [3] shows the voltage across a sample as a func- 
tion of time. One can see that the mean voltages obtains 
by time averaging are rather similar; the floating-kernel 
gives by 6% higher value than the TDGL theory. As 
one can expect from the transient nature of the floating- 
kernel correction, a more pronounced difference appears 
in the oscillation corresponding to passage of individual 
kinematic vortices. The floating-kernel approximation as 
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compared to the TDGL theory yields an increase in fre- 
quency of 17%. 

The time dependence in figure [3] was simulated un- 
der conditions corresponding to figure [T] In figure [1] one 
finds that the number of kinematic vortices in the TDGL 
approximation is six while in the floating-kernel approxi- 
mation there are seven. Since 7/6 — 1 ~ 17% we conclude 
that vortices move with similar velocities and the differ- 
ent frequency of oscillations follows from higher density 
of kinematic vortices. 

One can understand why vortex velocities are equal for 
both approximations in the case of an isolated Abrikosov 
vortex in the following way. The vortex velocity vy is 
given by the driving current density 3s and the friction 
coefficient ij as a. balance of the Lorentz and friction force 
Jg X z = $o^vy, where z is a direction vector either of 
the true magnetic field B or of the effective field Boff. As 
long as both fields are parallel, the vortex moves in both 
approximations with the same velocity. 

Although the velocity of kinematic vortices does not 
obey a simple law of balance of forces, applicable to iso- 
lated vortices, the numerical result suggests that their ve- 
locities are also only negligibly influenced by the floating- 
kernel corrections. During the time evolution of the sim- 
ulation, vortex core positions were traced and their ve- 
locities calculated. It was found that isolated vortices 
typically move about ten times slower than kinematic 
vortices. Since velocities of kinematic vortices change 
strongly with vortex positions^ it was not possible to 
compare the effect on velicities of the two approxima- 
tions with high accuracy. 



APPENDIX: FREQUENCY SUM RULE 

The conductivity is required to satisfy the frequency 
sum ruleiSii^ 



2 /■ , , , , ne^ 

— / dujiiea((jj,K) — , 

TT J m 



(24) 



where e, m and n are charge, mass and density of elec- 
trons. To satisfy this sum rule, a modification of the 
TDGL theory in the spirit of the two-fluid theory is neces- 
sary. Here we show that this two-fluid correction implies 
the floating-kernel approximation discussed in section HIl 



1. Sum rule violation in the TDGL theory 

First we show that the standard TDGL theory leads 
to a conductivity which violates the sum rule (P^ . In 
the standard TDGL theory one writes the total current 
as a sum of the supercurrent and the normal current, 



'TDGL — JGL -r- JN- 



(25) 



Neglecting the Hall effect, both currents are parallel to 
the electric field and the conductivity is a scalar given by 
the ratio, a — J/E'. It thus has two corresponding parts 

CTTDGL — CTQL + (TN- (26) 

The sum rule (1241) is satisfied in the normal state 
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- / dwRecrjv(w,k) = . (27) 



V. CONCLUSIONS 

Within the floating-kernel approximation we have dis- 
cussed the effect of the normal current on the motion 
of the condensate. It was shown that this effect is of 
transient nature and contributes only in time dependent 
systems with moving vortices. The most pronounced 
effects were found for kinematic vortices, which move 
much faster than isolated Abrikosov vortices. Particu- 
larly strong corrections were found near sample edges, 
where kinematic vortices are accelerated, increasing their 
speed nearly by an order of magnitude. 
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The superconducting component of mean Cooperon den- 
sity ns, mass m* and charge e* has an analogous sum 
over frequencies 

- / dwRe(TGL(w,k) = ''1^, (28) 
TT J m* 



The total sum rule (|24|) for cttdgl is thus violated. 



2. Two-fluid correction 

Assuming that formation of the condensate depletes 
the density of normal electrons, ^ n—2ns, the normal 
conductivity ought to be correspondingly lowered. 



= o-GL 
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(29) 



The sum over frequencies in the left hand side of ([M)) is 
then 



— / dwRe cr(a;, k) = ^ 
TT J m* 





" (30) 
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A sum rule similar to (l30l) was discussed in greater 
detail for the Meissner state, where a part of the weight 
due to superconducting electrons is covered by a singular 
term in the form of a Dirac S functioniiSr— Avoiding the 
6 function one arrives at the sum rulei^ 



2 f 

— / dwRe 17(0;, k) 

TT J 
+0 



^-^) ^. (31) 
n I m 



Here we assume similar structure for the mixed state. 
The only difference is that in the presence of vortices the 
conductivity ctgl of superconducting electrons is finite 
even at zero frequency. Its frequency dependence is not 
a 5 function but has a Drude form. 

In the pure limit, m* = 2m, the sum rule ^(IM is sat- 
isfied. As one can see, the sum rule (|24l) is violated in 
the dirty limit when m* 7^ 2m. This is in consequence 
of limitations of the theory used to derive the floating- 
kernel approximation. The derivation of 0] is based on 
the Kadanoff-Baym ansatz with the spectral function ap- 
proximated by the 5 function, therefore the renormaliza- 
tion of the Cooper-pair mass due to finite mean free path 
is not included in this approach. Briefiy, the floating- 
kernel approximation is justified only in the pure limit. 



3. Interaction of normal current with condensate 

We note that the floating-kernel approximation dis- 
cussed in section |IT] leads to the conductivity ([2^ . with 
current Here we approach the problem in the oppo- 
site way: starting from the conductivity (|29[) we arrive 
at the floating-kernel approximation. 

Let us require that the set of TDGL equations follows 
from the effective free energy-^ J- through 



r ( ^ - ie*cp ] ^ 
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Since our focus is on the spatial gradients, we neglect 
the non-linear relaxation of Kramer and Watts- Tobin. 
Indeed, the relaxation in the left hand side of (15^ results 
from ([2]), sending C — 0. 

In equations p2ll33p . the GL function t/i is normalized 
to the Cooperon density as ng — j^/'p and the sum rule 
uses the value averaged over space, ns = (IV'P)- We 
assume that the free energy has the superconducting and 
the normal parts, T = J-jv + ^s-, where the normal part 
is the same as in the normal state, therefore 
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In the TDGL theory the supercurrent reads 

Jgl = —Re ^ ( -ih\/ - — A I V- 
m* \ c / 



(35) 



(36) 



Since the vector potential A appears exclusively via the 
covariant gradient in the bracket, the current implies that 
in the TDGL free energy the kinetic energy takes the 
familiar form (l/2m*) |(-iW - (e7c)A) V-f- 

We have seen that the current ([55]) with the normal 
current added violates the frequency sum rule. Now we 
derive the kinetic energy assuming that current includes 
the two-fluid correction. According to the two-fluid con- 
ductivity (l29t. the total current reads 



^glE' + ( 1 - ^ j ajvE' 



'GL 



2ns 



n 



e - / e m \ 

= — Re ip { -ihV - —A Jty MA + Jw- (37) 

m* \ c en J 

We have included the correction term in the supercur- 
rent, because it is proportional to the condensate density. 
According to ([55)) and ([SSt it thus cannot result from the 
variation of the normal free energy J^^v- 

From equations (|33l) . (|34|) and (|37l) one flnds 



SA 



Re V -iW 



en 



3n i'- (38) 



Integrating relation ([38)) over the vector potential one 
finds the superconducting free energy of form 



^5 



1 



2m* 



-thV A J 



N 



+ a|VP + -/3|V'|'. 



(39) 



The GL function thus enters the superconducting part 



Of course, the integration provides only the kinetic en- 
ergy which has to be rearranged with integration by parts 
into the square of covariant gradients. The terms inde- 
pendent of A represent initial conditions for the integral 
and are taken from the standard GL theory. 

With free energy (15^ , equation ([5^ is identical to the 
fioating-kernel approximation ([2]) in the C — ?• limit. 
It should be noted that derivation of equation ^ from 
microscopic theory was also carried to terms linear in 
Jjv, therefore additional terms quadratic in the normal 
current might appear. 

To summarize this appendix, we have shown that the 
TDGL theory violates the longitudinal /-sum rule. To 
restore the sum rule one must reduce the normal cur- 
rent which corresponds to the interaction term between 
the normal current and the condensate. In this way 
one recovers the floating kernel-approximation from phe- 
nomenological arguments. 
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